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Abstract 


The paper aims at estimating and forecasting international 
tourist arrivals to Cambodia during the time interval of 
2000m1 to 2017m7, covering 209 of monthly 
observations. To find out factors affecting tourist arrivals, 
simple OLS and 2SLS with instrument variable regression 
are applied, on the one hand. On the other hand, several 
time series models of ARIMA (p, d, q), GARCH (s, r) and 
the hybrid of ARIMA(p, d, q)-GARCH(s, r) are employed 
to forecast tourist arrivals in line with AIC and BIC in 
selecting the best modified models. The empirical results 
primarily reveal that tourist arrivals are affected by 
exogenous factor, say exchange rate, dummy factors such 
as the AEC, global finical crisis, national election and 
Cambodia’s e-Visa. With regard to forecasting stage, the 
result indicates that tourist arrivals are shocked by time 
trend in the past period, say time (t-1). The trend is 
furthermore reduced due to the time lags, say time (t-2, t- 
3) as shown in the parameter coefficients of AR. GARCH 
(1, 1) model suggests that the short run persistence of 
shocks lies in the gap of 0.04 whereas the long run 
persistence lies in the gap of 0.94. Additionally, AIC and 
BIC propose that ARIMA(3, 1, 4) and the hybrid of 
ARIMA(3, 1, 4)-GARCH (1, 1) are the best model to 
predict the future value of tourist arrivals. The RMSE and 
U index obtained from measurement predictive accuracy 
reveal that long run 1-step ahead forecasting of 2013m12 
to 2017m7 is produced the smallest predictive error 
amongst the others. Thus, it has more predictive power to 
apply long term ex-ante forecasting. 
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Forecasting, ARIMA (p, d, q)- GARCH (s, r) Model, 
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1. Technical Observation 


Prediction or forecasting is generally considered as an art of anticipation or estimation 
any future event and/or value. In the context of economic and financial time series 
analysis, it somehow takes into account the prediction methods due to econometric 
models in line with statistical inferences. It is knowingly separated by two main 
categories that so-called in sample and out of sample forecasting or say ex-post and ex- 
ante forecasting. Good performance in out-of-sample prediction is viewed as the acid 
test for a good forecast model (Kunst 2012). It reflected the facts of any econometric 
models which is perfectly and methodologically adopted. In this case, diagnostic tests 
are employed conventionally. Forecasting in macroeconomic or financial data is widely 
acknowledged since it has played an important catalyst for policy makers as well as 
financial trader to set up the policy in achieving growth, development of the country and 
to gain profit from market speculation respectively. Prediction in general vastly meet 
the maturity. Many different approaches, both linearity and nonlinearity, due to the 
combination method of mathematics and/or statistical inferences are applied, (Wang 
2016), (Sjo 2011), (Chia-Lin Chang 2009) and (Elzbieta F. and Mirostaw Ga 2004). 
Time series models can be by definition giving a reasonable benchmark to evaluate the 
ranging value of forecasting based on periodical step ahead and/or full and/or sub 
sample observation relevant to the pure explanatory power of historical behavior of the 
series if the methodological assumption is detected and not violent the assumptions 
within the models such as the Box-Jenkins methodology of ARIMA (p, d, q), (Box and 
Pierce 1977). 


In financial and economic time series estimation and prediction, the most common 
models which were typically and frequently employed are autoregressive conditional 
heteroskedasticity or so-called the ARCH model, (Baum 2015). Using an ARMA 
processes with up to two lags and variance with one of GARCH, EGARCH or TARCH 
processes with up to two lags, (Jansky 2011) evaluated several hundred one-day-ahead 
of VaR forecasting models. GARCH process as the best conditional volatility process 
for the analyzed time series, stated by the above author. It helped improvements, 
including a no prior assumption on the distribution of the log returns, which proved to 
be a step in the right direction. Consequently, time-variation in volatility 
(heteroskedasticity) is a common feature of financial data. The most straightforward 
way to measure the heteroskedasticity is to estimate the time-series of variances on 
“rolling samples”, (Chen 2013). This can model by considering heteroskedasticity. Yet, 
in the context of tourist arrival forecasting, it can be defined that as time-varying 
conditional variance has both the AR and MA components, it leads to the generalized 
ARCH (p,q) (GARCH (p,q)) of Bollerslev (1986), (Chia-Lin Chang 2009). In order to 
manage international tourism growth, therefore it is essential to model sufficiently 
tourist arrivals and their association due to volatility and autoregressive model and 
forecasting in the study of (Chia-Lin Chang 2009). Mamula (2015) suggested that 
although the diagnostics for the selected models reveals that the proposed models do not 
significantly differ, it can be concluded that the multiple regression model performs a 
highly accurate forecasting of tourist arrivals. 


Tourism demand forecasting currently and widely employed the variety of forecasting 
methods, running from the most rudimentary approaches to the more complex, (Ramos 
2014). Tourist arrivals and expenditure (receipt), in both aggregate and per capita forms, 
are commonly used to measure tourism demand in empirical research, (Song et al. 
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2010). In line to its modelling and forecasting, most of the empirical studies have 
focused on conventional approaches of forecasting performance toward measurement 
predictive accuracy, those models are included but not limited univariate ARMA and 
ARIMA based models, GARCH or the hybrid ARMA-GARCH, ARIMA-GARCH, 
seasonality components as well as nonlinear models, (Chia-Lin Chang 2009), (K.-Y. 
Chen 2011), (Andrew Saayman 2015), (Robert R. Andrawis 2011), (Chu 2009), (Shan 
2002) and (Witt 2003). Tourism demand and volatility modelling and forecasting 
(Suhejla Hoti n.d.), (Chikobvu 2017), (Louis 2015), (McAleer 2005) and (T. K. Chia- 
Lin Chang 2011). 


With respect to few exogenous and dummy factors, (George Agiomirgianakis 2014), 
found the negative relationship between exchange rate volatility and tourist inflows into 
Turkey. The study of exchange rate volatility and tourism demand, (Webber 2001), 
(Yap 2012) and (Chang 2009). In contrast, (Crouch Geoffrey I. 1993) underlined that 
exchange rate volatility is a contributing factor to tourist arrivals. Both the moving 
average and the high and low measures of volatility have proven to have a significant 
effect to tourist arrivals. In addition, the impact of financial crises on tourism demand is 
less significant. Ensuring the safety and health of tourists is the key to maintain demand 
for inbound tourism, (Yu-ShanWang 2009). Since the 2008 global financial crisis and 
resulting recession, many countries have been following unconventional monetary 
policies. Some findings highlight the importance of government economic policy in 
stimulating international tourism demand through its impact on the economy, (Jewoo 
Kim 2016). 


The central purposes of the study are firstly to find out the factors which might affect to 
tourist arrivals in Cambodia. It aims secondly at modeling and forecasting tourist 
arrivals through the autoregressive and volatility approach. Within this stage, the study 
accordingly adopts 1-step ahead of out of sample forecasting to check the measurement 
predictive accuracy as well. Indeed, the remainder of the paper is organized as follows. 
The 1* session is to present the facts of time series prediction and factors affecting to 
tourist arrivals whereas the 2" one is to design the methodology and data calculation. 
The 3 and last session is to interpret the empirical outcomes and concluding remarks 
along with suggestions. 


2. Methodology and Data Calculation 


2.1. Data Collection and Calculation 


Responding to our central objectives, international tourist arrivals in Cambodia, simply 
noted as (TA) is employed. Monthly observations of TA from 2000m1 to 2017m7 are 
extracted from ministry of tourism, (MoT, 2017). Exchange rate (EX) is imported from 
the CEIC data manager. The study transferred the sample observations of TA into the 
nature of log return using the formula as follows: y; = In(TA,) — InCTA__1) for the 
forecasting stage. Therefore, the descriptive statistics of both TA and EX with and 
without taking logarithm function are demonstrated in table | as bellows: 


Table 1 shows the different types of descriptive statistics such as mean, minimum and 
maximum as well as LM test so on. With the total observations of 209, Shapiro-walk 
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statistics' and normality test show that the series TA and EX are not come from normal 
distribution assumption. LM test for ARCH effect somehow reveals the rejection of the 
null hypothesis at 1 % level of significance. It statistically defines that the series are 
contained ARCH effect due to the null hypothesis of no ARCH effect”. Skewness and 
Kurtosis statistics give insights into the shape of the normal population distribution. 
Skewness essentially measures the relative size of the two tails or say the normal 
distribution has a Skewness of 0 whereas Kurtosis is a measure of the combined size of 
the two tails. Its value often compares to kurtosis which is equal to 3 or greater. As the 
result, both statistics indicate that Skewness value is ranged in 0.64 and -0.39 and -0.51 
and -0.59 respectively whereas Kurtosis is 2.5 and 2.06 for nature and logarithm data of 
TA, respectively. It statistically indicates that the series, TA; and In(TA,) is non-normal 
distribution as similar as denoted in the Shapiro-walk test as well. 


Table 1: Descriptive statistics of tourist arrivals to Cambodia, TA 


Description TA, In(TA;) EX; In(EX;) 
Observations 209 209 209 209 
Percentiles (50%) 173112 12.0617 4104.482 8.3198 
Mean 204768.4 11.9678 4089.954 8.3159 
Standard Deviation 135157.1 0.7782 115.868 .02856 
Min 30485 10.3249 3811.999 8.246 
Max 611534 13.324 4305.74 8.368 
Variance 1.83e+10 0.6056 13425.46 0.0008 
Skewness 0.6354 -0.3937 -0.5136 -0.588 
Kurtosis 2.5012 2.0681 3.004 3.088 

5.4 4.65% 4.22% 4 52 

prep Wallk Test (0.0000) (0.0000) (0.0000) (0.0000) 
140.01 *** = 174,29%% 181.79%** 183.19%*# 

EN UteS tora CED Ciea (0.0000) (0.0000) (0.0000) (0.0000) 


Source: Author’s estimates 

Note: The sign notification of *, ** and *** refereed to the statistical significance of 10%, 5% and 1%. LM test for 
ARCH considered 1 degree of freedom in order to test after a regression of its own trend. The value inside the 
parenthesis is the p-value. 


2.2. Estimation and Prediction Approach 
2.2.1. Baseline Regression Model of Tourist Arrivals in Cambodia 


To estimate the factors affecting tourist arrivals in Cambodia during the period 
observations, let’s consider a sample regression equation in line with time trend effect, 
(TT,) as follows: 


TA, =C +y[X':] + t|TT;| + u,,07.~N(0, 1) (1) 


Where TA; is an explained variable and it is denoted as international tourist arrivals. u; 
is an error term and the constant term, (c). X’; is a matrix set of explanatory variables. 
The study employs exchange rate as the main explanatory variable. It is known partially 
as tourism price effect. TT; is a matrix set of time trend effect during the observed 
periods. To overcome the time trend effect on the model, dummy (binary option) 
variables which take into account number | for the determined period and 0 otherwise, 


' Shapiro - Wilk and Shapiro - Francia tests for normality, https://www.stata.com/manuals13/rswilk.pdf 
> Arch - Autoregressive conditional heteroskedasticity (ARCH) family of estimators, STATA Journal, 
http://www.stata.com/manuals 13/tsarch.pdf 
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will be adopted. Those dummy variables are included the national election in 2003, 
2008 and 2013, the global financial crisis in 2008 and 2009, milestone of the ASEAN 
Economic Community (AEC) in 2015 and Cambodia first e-Visa launching from 2006 
to present. These variables believe to have the strong impact to travel decision. Yet, the 
study applies OLS estimator with robust standard error (SE) and 2SLS with instrument 
variable (IV) to estimate the equation (1). 


Furthermore, to model and forecast the volatility and the time trend effect of tourist 
arrivals individually, the study applies ARIMA (p, d, q), GARCH (s, r) and the hybrid 
ARIMA (p, d, q)-GARCH (s, r) model. Therefore, ARIMA (p, d, q) is modeled as 
follows: 


2.2.2. ARIMA and ARIMA-GARCH Model 


Most of time series data are econometrically affected by either autoregressive process 
(AR) or moving average process (MA). Let’s consider, in one part, an ARIMA (p, d, q) 
with parameter order of p and q in line with d (order of integrated, D) as follows: 


Ve = A + OVp-1 + O1Vp-2 ++ PyVi-p + OyUp_1 + PyUp-1 Fo + Pglt-q (2) 


It is noteworthy that equation (2) can be derived with different data and lag operators or 
AR (p) and MA (q) process. Hence, we get: 


Ve = A + yee OyVt—p + eT Pqlt-q + TT (3) 


Or using the backshift notation with By, = y,_,, the above equation, (3) can write as 
follows: 


(1-0,B-+— OB? )y, = (1-@iB—-— BDU 4) 


Furthermore, using the first different of series, W, = yz — yz_1 or W; = (1 — B) yz, the 
specific general form of ARIMA (p, d, q) is equated as follows: 


Dy(B)(1 — B)*y¥¢ = Gq(B) ur (5) 


Determining the parameter in different equation is a must in ARIMA (p, d, q) due to 
Box-Jenkins methodology, (Box 1970) and (Box and Pierce 1977). Accordingly, 
checking stationary or unit roots is essential and important to decide the parameters of 
all elements in the model. The study employs two popular methods of unit roots test, 
namely anaugmented Dickey—Fuller (ADF) tests, (Dickey 1981) and Phillips- 
Perron (PP) test, (Phillips 1988). The study tests these two tests with and without trend 
and intercept. It is noted that these tests contain the null hypothesis of having unit root 
or meaning that the series is non-stationary. As the result, the empirical outcomes 
demonstrate in table 2. Without trend and intercept, both series are non-stationary at 
level and they are stationary at first different, say I(1) due to the converting. Conversely, 
with trend and intercept, the series is stationary at level, say I(0). As the result, to 
control the interaction of log likelihood not to be concave iteration in post estimation, 
logarithm series with first difference will be employed. Next, the study presents a brief 
description of GARCH (s, r) model. 
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Table 2: Unit roots test analysis of TA and EX 


= TAL In(TA,) EX, In(EX;) 
Desenpton. SN IT NT IT NT IT NT IT 
At level, I(0) 
ADF -1.76 -5.32 -1.82 -4.88 -2.73 -2.64 -2.74 -2.63 
(0.4003) (0.0001) (0.3722) (0.0003) (0.0686) (0.2631) (0.0682) (0.2648) 
PP (Z(rho)) -6.28 -59.86 -5.01 -53.24 -14.69 -17.16 -14.51 -17.01 
(0.3900) (0.0000) (0.3763) ~— (0.0001) 0.0354 0.1251 0.0358 0.1271 
At first difference, I(1) 
ADF -12.53 -12.5 -11.98 -11.96 -11.17 -11.17 -11.15 -11.15 
(0.0000) — (0.0000) (0.0000) (0.0000) ~— (0.0000) ~— (0.0000) ~—_ (0.0000) ~— (0.0000) 
PP (Z(tho)) -160.48 -160.4 -146.06 -146.05 -147.75 -147.86 -147.41 -147.53 


(0.0000) (0.0000) (0.0000) (0.0000) (0.0000) (0.0000) (0.0000) (0.0000) 


Source: Author’s estimates 

Note: The statistical value in the parenthesis is p-value. P-value of 0.0000 indicated the statistical significance of 
rejection the null hypothesis at 1%. NT is referred to the estimation without trend whereas IT with trend and 
intercept. 

With regard to volatility approach, ARCH model introduced by (Engle 1982) and 
generalized ARCH, the so-called GARCH (Generalized ARCH) by (Bollershev 1986) is 
used to investigate the volatility effect in the series, both low and high frequency data. 
The models widely adopt in various branches of econometrics, particularly in financial 
time series analysis. To estimate and forecast TA, a standard GARCH (1, 1) with no 
regressors in the mean and variance equations is proposed. Therefore, the model is 
equated as follows: 


Y; = OX", + Et, Et = Of¢Z¢ (6) 
O77, =wt ae. 1 + Bory (7) 


Since 07; is the one-period ahead forecast variance based on past information, it is 
called the conditional variance. €7,_, and o7;_; are an ARCH and GARCH term 
respectively. Parameter testing states that z, is standardized residual returns (i.e. iid 
random variable with zero mean and variance. For GARCH (1, 1), the constraints 
a, B = 0 and w > 0 is needed to ensure that a”; is strictly positive, (Suliman 2011). As 
the result, from equation (5) and (7) we can derive the hybrid of ARIMA (p, d, q)- 
GARCH (1, 1) model as follows: 


0,(B)(1— B)*y, = Pq(B)ur + +a€*;_4 + Bo ¢-4 (8) 
2.2.3. Measurement Predictive Accuracy 


Since the paper aim at forecasting tourist arrival from post estimation and out of sample 
prediction, the measurement predictive accuracy is adopted. Hereafter, let’s assume 
forecast sample is j = T + 1,T + 2,...,T +h, and denote the actual and forecast value 
in period t as y; and ¥; , respectively. The study uses two types of error predictive 
method, namely root mean square error (RMSE) and Theil’s inequality index (U) to 
measure. Furthermore, it is used to specify the best model for the purpose of long run 
ex-ante prediction. Therefore, the forecast evaluation measurements, RMSE and U 
define as follows: 


RMSE = [X72 (9 —y,)2/h (9) 
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[Bet es@e-v0?/n 


Theil’s U = (10) 


[amit oo2/n+ JETth, ow?/n 


Shortly, the study employs ARIMA (p, d, q), GARCH(s, r) and the hybrid of ARIMA 
(p, d, q)-GARCH (s, r) to model and forecast tourist arrivals in Cambodia from monthly 
observation of 2000m1 to 2007m7. The RMSE and U are employed to measure the 
predictive accuracy from out of sample forecasting. 


3. Empirical Result and Discussing 


This session aims at interpreting the empirical outcomes from post-estimation of the 
baseline regression equation, equation (1) and forecasting models from ARIMA (p, d, 
q), GARCH (s, r) and the hybrid of ARIMA (p , d, q)-GARCH (s, r). Indeed, the study 
presents some calibration of tourism trend and development since 1993 to present 
toward numerical sources and graphical illusions. 


3.1. Calibration of Cambodia’s Tourism Industry 


Table 3 shows the tourism trend in Cambodia since 2013 till February 2017 and its 
growth rate year-on-year. From 2016 to the first quarter of 2017 (Q1), tourism growth 
rate approximates 12.1% comparing to those of 2015 and 2016 which presents 5% and 
6.1% for 2015/2014. It demonstrates from year to year that tourism trend is increased 
considerably due to on the one hand government considers tourism industry as one of 
the central sector in contributing to growth and development. On the other hand, 
Tourism Development Strategic Plan 2012-2020 is implemented with the goal of 
attracting tourist throughout connectivity, safety and security, marketing and facilitation 
of tourist transportation, (MoT 2012) for example. In 2013, international tourist arrivals 
account over 4 million people while it reached up to almost 5 million people in 2016. It 
reflected the gap of slowing growth, as it presents just 1 million people coming to visit 
Cambodia for almost 4 years. According to the research, only 80% have visited 
Cambodia one time and 20% is more than one. Furthermore, lack of infrastructure, 
especially the sewage facilities, inadequate accommodations and facilities, personal 
security issues, and accessibility to secondary destinations are items which need 
immediate attention, (Paul Leung 2017). 


A decade afterward of UNTAC intervention, international tourists reached up to around 
700 thousand visitors in 2003, it rises almost 500% comparing to 1993. Currently, the 
trend is skewed slightly for more than 5 million visitors which equaled to one three of 
the Cambodian population. For instance, the industry is knowingly shocked by both 
Asian crisis in 1997, the growth rate has diminished by 16% and in global crisis in 
2008, has reduced by around 11%, (Figure 2). Tourist receipts are augmented 
dramatically, accounted 3.2 million USD in 2016 as of that in 1995 attracted 1 million 
USD. Most of inbound tourists coming to Cambodia spent 6 to 7 days for staying with 
the pocket payment of 640$ per tourist (data calculated from the tourism statistics 
report, April 2017). 
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Table 3: International tourist arrivals to Cambodia, 2003 - 2017 


Change (%) 


15/— 16/ 17*/ 


Months 2013 2014 2015 2016 2017 14 15 16 
Ql 1,172,072 = 1,267,922 1,307,836 1,342,477 1,025,521 3.1 2.6 

January 404, 106 442,045 460,577 466,086 532,206 4.2 1.2 14.2 
February 385,760 425,801 430,207 448,468 493,315 1 4.2 10 
March 382,206 400,076 417,052 427,923 18.3 2.6 

Q2 920,527 933,446 994,154 1,018,455 0 65 2.4 

April 327,000 332,690 361,139 367,684 17.9 1.8 

May 292,115 300,302 314,748 320,601 253 1.9 

June 301,412 300,454 318,267 330,170 20.1 3.7 

Q3 964,612 998,690 1,044,880 1,147,483 0 46 98 

July 338,761 340,091 364,325 395,761 19.2 8.6 

August 342,064 347,211 366,096 406,214 16.4 11 

September 283,787 311,388 314,459 345,508 16.9 9.9 

04 1,152,954 1,302,717 = 1,428,361 — 1,503,297 0 96 5.2 

October 334,410 390,637 408,922 414,077 14.9 1.3 

November 386,737 411,501 444,640 477,686 16 7.4 

December 431,807 500,579 574,799 611,534 12.9 6.4 

Total 4,210,165 4,502,775 4,775,231 5,011,712 = 1,025,521 6.1 5 12.1 


Source: Tourism Statistics Report, Ministry of Tourism (2017) 


The strong commitment and policy base development of ministry of tourism and the 
related institutional state is putting in place for strong support. Well administrative and 
foreign policy toward China has brought Cambodia one of the most attractive place for 
China tourists. From 2010 to 2017, there are 7 million China tourists travelling to 
Cambodia. From which Cambodia is a home of both natural and cultural heritage, there 
are places to enjoy and leisure. Currently, the efforts have been made to attract 
additional arrivals by establishing more direct flights and introducing new initiatives 
such as the “China Ready” initiative and joint tour packages. 


Figure 1: International tourist arrivals from 1993 to 2016 


6.000.000 60,0 
49,4 50,0 
5.000.000 
40,0 
4.000.000 
30,0 
24,4 
3.000.000 20,0 
10,0 
2.000.000 
0,0 
1.000.000 ; | 10,0 
Qo =F | | i | | | I -20.0 
ntrnonw~anaeorantuve tr wnanonnantne : 
DanagaananngonooermrceocmcmcUOCUCUCMOCUCUCcUOCUCUOUOUCUCUW CUP BR hOB OD OD ae se SS 
DARDADAARARAASCSSESESSESESCSESSSCOSCSSESSS 
oD so Ds Do Ds DC os BS Co DC 
ma Number of tourist arrivals —— Change (%) 


Source: Tourism Statistics Report, Ministry of Tourism (2017) 
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Figure 2: Average length of stays and international tourism receipts, (1995 — 2016) 
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Figure 3: Share of monthly tourist arrivals from 2013 to 2016 
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Source: Ministry of Tourism, Cambodia (MOT), (2017) 


3.2. Estimation Outcomes from the Baseline Regression Equation 


To estimate factors affecting tourist arrivals in Cambodia, the study employs monthly 
observations of international tourist arrivals as an explained variable, exchange rate as 
the main explanatory variable and some dummy variables such as the global financial 
crisis, the national election, AEC and Cambodia’s e-Visa as an exogenous variables. 
The empirical results show in table 4. With 209 of sample observations, all opposed 
models indicate the statistical significance with 1% level in line with F-statistic value. 
Thus, they are perfectly and correctly modified. From model (1) to (4), the study tests 
baseline regression through OLS with robustness SE whereas model (5) to (6) apply 
2SLS with IV regression. Model (1) and (2) estimate the exchange rate only in the 
before-after the global financial crisis without controlling dummy variables. From 
model (3) to (6), the study uses exchange rate in the whole sample to estimate with and 
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without dummy variables. In addition, model (5) imports all dummy variables as the 
instrumental variables and adopts 2SLS instrument variable regression. 


It is reflected that exchange rate before the global financial crisis are negatively affected 
to TA with the statistical significance of 1%. This trend is showdown tourist arrivals 
resulted of rising exchange rate. Conversely, exchange rate after crisis and exchange 
rate in the whole sample, model (2) to model (6), are positively associated with TA. 
This is likely due to the facts that Cambodia has pegged its own currency to US dollar 
(USD), its appreciation and depreciation might not reflect strongly to tourist decision. 
Or meaning that Cambodia’s currency is allowed the value to be fluctuated due to 
supply and demand of market. The positive is given the facts of effecting and 
supporting the importance of maintaining a relatively stable exchange rate to attract 
international tourist arrivals. In addition, tourist arrivals to Cambodia is not sensitive to 
currency shock since Cambodia is able to maintain the appreciation and depreciation in 
the relative gap. 


Table 4: Factor effecting tourist arrivals in Cambodia 


Baseline Regression Equation 


Tourist arrivals 


Model (1) Model (2) Model(3) Model (4) Model(5) Model (6) 
Pechanceats 14.526*** = 4,758*** 36.16*** 4.758 *** 
(10.46) (4.69) (9.18) (4.63) 
Exchange rate x = 27 *** 
before-crisis (-19.67) 
Exchange rate x Li25*¢* 
after-crisis (20.29) 
Dummy variables 
National election ree 108) 
3 Miag 1.148*** 1.148*** 
(16.40) (17.70) 
Global financial -0.392*** -0.392*** 
crisis (-5.96) (-4.57) 
0.569*** 0.569*** 
nue (9.64) (7.38) 
Consiantieen 12.55*** = 11.43*** — -108.83*** = -28.39*** — 288.8 ** -28.39** 
(363.23) (216.00) (-9.43) (-3.37) (-8.81) (-3.33) 
Tumnber ot 209 208 209 209 209 209 
observations 
ee 386.77*** ALL.77*** — 109.51*** = 163.16*** = 84.25% ** 163.16*** 
(0.0000) (0.0000) (0.0000) (0.0000) (0.0000) (0.0000) 
Adjust R? n/a n/a 0.2843 0.8007 n/a 0.7958 


Source: Author’s estimates 

Note: Robust t statistics in parentheses and * p<0.05, ** p<0.01, *** p<0.001. OLS is estimated in line with robust 
standard error (SE). IV is estimated through instrument variable (IV) technique or so-called Instrumental variables 
(2SLS) regression. 2SLS (1) used all dummy variables to be an instrument variables. n/a defined not available 
information. 


With respect to dummy variables, it is revealed the empirical outcomes as expected. The 
global financial crisis is negatively linked to tourist arrivals. It is likely occurrence to 
some empirical studies that a significant slowdown in the Turkish foreign active tourism 
during the global crisis (Kudret Gul 2014). In addition, (José F. P-R. et al. 2016) found 
that the proposed model is appropriate for explaining the changes in the market 
positions caused by the economic crises. According to (UNWTO 2013) stressed that the 
2008-2009 global economic crisis severely affected international tourism, causing in 
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2009 a decline of 4% in international tourist arrivals and a decrease of 6% in 
international tourism receipts in 2009. The crisis actually caused the first serious 
downturn faced by international tourism in decades, a sector accustomed to a long-term 
average growth rate of about 4% a year. The World Bank also stated that the financial 
crisis has cut access to loans in advanced and developing countries, pulling investment 
out of poorer nations and reducing consumer spending. Hence, reducing consumption 
and investment, showdown tourist travel, then. 


Furthermore, Cambodia’s e-Visa and the AEC are found to be positively related to TA. 
Launching and adopting such the technology innovation brought Cambodia as an easy 
place to access to travel, particularly applying for short term Visa. This e-Visa seems to 
have been embraced. Certainly, it has reduced the time and expense required in securing 
official permission to travel to the sub-continent. Hence, Cambodia is one of the easiest 
countries in the world to emigrate to, visa-wise. The positive of the AEC is undoubtedly 
since AEC will bring not only trade and investment flow but service, labors as well as 
tourist arrivals with visa exemption for some countries within the nations. The study 
furthermore stressed that the national election is also positively associated to TA but it 
is shown a small proportion, meaning that it is not existed the strong impact. It is 
opposite to (Ghana | Bennet Otoo 2016) harassed that the word ‘’elections’’ is often 
surrounded by a general stigma of fear, chaos, and anxiety. In every part of the world, 
the electioneering period is a brief period of a dip in almost every sector of life. A lot of 
activities are put on hold and investors/businessmen are reluctant to travel or do 
business in such countries at this time. 


3.3. Estimation Outcomes from ARIMA and ARIMA-GARCH Models 


To estimate and forecast tourist arrivals in line with TA variable individually, the study 
employs TA as the nature of the logarithm with the first different, I(1) to overcome the 
stationary process. Baseline regression between TA and TA at lag 1 is adopted to 
capture the pattern of TA at level, I(0) and first different data, I(1). The result shows in 
figure 4. It somehow defines that the simple regression with time trend effect and scatter 
plotting of sample data draw in the gap of 10 to 14 as converting to logarithm function 
and in the gap of -4 and 4 due to taking the first different data. TA at time (t-1) is 
positively and significantly associated to TA at time (t). It reflects that the more tourist 
arrivals in the past year, the more they visit in the present year. This gives the idea that 
tourists return to the country of visit due to level of pleasant and satisfying. According 
to research, there is 20% of tourist arrivals in Cambodia have been visited more than 
one time. 


To apply Box-Jenkins methodology, detecting the random order and stationary process 
is though applied in the previous session, the study employs autocorrelation (AC) to 
determine the best fit parameter of AR and partial correlation (PAC) for MA parameter. 
As the result, it is shown in figure 5 as bellows. 
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Figure 4: Post-estimation in line with level and log differential data 


Source: Author’s estimates 


Figure 5 specifies that ACF contains significant values in the first three lags, while the 
PACE exhibits decay in the form of an approximate damped. As the result, parameter of 
AR could be tested by 1, 2, 3 and MA is tested from 1, 2, 3 and up to 4. It suggested as 
an appropriate specification. The study selects the parameters of AR and MA due to 
graphing AC and PAC. The result suggests to adopt its parameter differently. Therefore, 
its empirical results show in table 5. The study proposes a post estimation from 
univariate time series models, namely ARIMA, GARCH and hybrid of ARIMA- 
GARCH model with the interaction of time trend (t). It employs full sample of 2000m1 
to 2017m7, covering 208 observations. 


Figure 5: Autocorrelation (AC) and partial correlation (PAC) at log different data 


Source: Author’s estimates 
Table 5 shows that constant term (c) and time trend is still insignificant and keep the 


sign constantly. More interestingly, it is such a straight forward that constant term (c) is 
insignificant. It is not a must to drop the constant term in the models*. Most of AR and 


3 Usage Note 23136: Understanding an insignificant intercept and whether to remove it from the model, 
URL: http://support.sas.com/kb/23/136.html 
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MA coefficients are negatively and positively affected to TA. In model (1), say ARIMA 
(1, 1, 1) demonstrates that MA(1) is positively affected to TA whereas AR(1) reveals 
insignificant association. The positive of MA reflects the statistical method that tourist 
arrivals (TA) contains an MA term in the model. Coefficient of AR and MA of model 
(2), (3) and (4) show the significant relationship to TA. Coefficient of AR(1) is 
positively affected where AR(2) and AR(3) are negatively related to TA. The positive 
relationship of AR(1) to TA reflects the facts of tourist arrival at the present time, say 
time (t) is impacted by tourist at past time, say (t-1) but it is reduced by time destructive. 
It means that when time trend is reduced due to AR coefficient, tourist arrivals will 
reduce as it is shown a negative impact, AR(2) and AR(3). 


Table 5: Tourist volatility, estimation of full sample observations 


GARCH ARIMA (p, d, q)- 


Models eo ee (s, 1) GARCH (sr) 
() (2) (3) (4) (5) (6) (7) 
Constant term 0.06 0.07 0.05 0.06 0.05 0.06 0.05 
(c) (0.55) (0.46) (0.76) (0.59) (0.52) (0.58) (0.42) 
Sigma 0.15*** 0.15*** 0.15*** O.15*** 
(17.62) (16.97) (18.48) (17.12) 
AR parameters 
0, -0.2 0.24*** -0.22 
(-0.77) (3.30) (-0.82) 
@> -0.28*** 
(-3.70) 
@3 -0.13** -0.13** 
(-1.79) (-1.71) 
MA parameters 
Pi 0.45** 0.46** 
(1.93) (1.96) 
D4 0.25*** -0.08*** 0.13** 0.13** 
(3.64) (-1.00) ( 1.94) (1.94) 
ARCH & GARCH parameters 
a 0.04 0.0004 0.04 
(0.93) (0.00) (0.92) 
B 0.94*** 0.19 0.94*** 
(12.16) (0.00) (12.76) 
w 0.0002 0.0177 0.0004 
(0.19) (0.00) (0.49) 
Testing Parameter Coefficient in GARCH (r, s) 

a,p =0 672.84*** n/a 703.70*** 
a+p=0 481.62*** nla 511.12*** 
Time Trend Yes Yes Yes Yes Yes Yes Yes 

AIC -199.9 -200.92 -199.26 -193.47 -193.36 -189.52 -201.78 

BIC -183.22 -184.24 -182.57 -176.78 -176.67 -166.16 -178.22 

Wald chi2 21.58***  IOTLEEF 14.95% ** 6.49* 0.17 6.07* 18.10*** 
Observations 208 208 208 208 208 208 208 


Source: Author’s estimates 

Note: Model 1, 2, 3 and 4 refer to ARIMA (1, 1, 1), ARIMA (1, 1, 4), ARIMA (2, 1, 4) and ARIMA (3, 1, 4) 
respectively. Model 5, 6 and 7 denote GARCH (1, 1), ARIMA (3, 1, 4) - GARCH (1, 1) and ARIMA (1, 1, 1) — 
GARCH (1, 1) respectively. The sign notification of *, ** and *** refereed to the statistical significance of 10%, 5% 
and 1%. 


Therefore, the most suitable model to adopt an out of sample forecasting is model 4, say 
ARIMA (3, 1, 4) due to the lowest statistical value of AIC and BIC and significant level 
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of the Wald chi-square, which is approximated 6.49. Yet, ARIMA (3, 1, 4) could be 
considered to model with GARCH (1, 1) as the hybrid model, say ARIMA (3, 1, 4)- 
GARCH (1, 1). Consequently, ARIMA (1, 1, 1)-GARCH (1, 1) and model 7 takes into 
account. From model (5) and (6), the study applies GARCH (1, 1), ARIMA (3, 1, 4)- 
GARCH (1, 1) and ARIMA (1, 1, 1)-GARCH (1, 1). It shows that AR(3) remains 
negatively associated with TA where AR(1) is not significant impact. In line with 
ARCH and GARCH coefficient specify only # is significant impact. The insignificant 
of @ reveals that volatility shock today of tourist arrivals is not fed through into next 
period’s volatility. More importantly, GARCH(1, 1) model of tourist arrivals to 
Cambodia suggests that the short run persistence of shocks lies in the gap of 0.04 while 
the long run persistence lies in the gap of 0.94. As the second moment condition, a + 
B = Ois satisfied. 


3.4. Estimation Outcomes of Measurement Predictive Accuracy 


In this session, we aim at measuring the forecasting error from out of sample prediction 
due to the lowest value of RMSE and U index. This adopts the appropriated models 
throughout ARIMA (p, d, q) and GARCH (s, r) with one step (1-step) ahead obtained 
from table 5. The first 1-step ahead applies the post estimation from sample 
observations of 2000m1 to 2013m12. Forecasting out of sample is afterward considered 
from 2013m12 to 2017m7. The second 1-step ahead starts from 2014m12, the third 1- 
step ahead from 2015m12 and the fourth 1-step ahead from 2016m12 to 2017m7, 
respectively. As the result, the statistical value of RMSE and U report in table 6 and the 
comparison of the models illustrates in figure 6. It demonstrates that from the model (1) 
to model (5), RMSE and U index is beyond and insight into 1-step ahead as it produced 
the smallest error amongst others. It is likely due to the facts that the models are 
perfectly fitted with long period ahead rather than the nearest period. 


Figure 6: One-step prediction of residuals from the post estimation 


Ki 
inh 
7 AN | 


xb prediction, one-step 


Mn) | 


T 
2000m1 =2002m1 2004m1 2006m1 2008m1 2010m1 2012m1 2014m1 2016m1 2018m1 


ARIMA(1, 1), one-step 
ARIMA(2, 1, 4), one-step 
ARIMA(2, 1, 4), one-step 


ARIMA (1, 1, 1) - GARCH(1, 1), one-step 
ARIMA (3, 1, 4) - GARCH(1, 1), , one-step 


Source: Author’s estimates 
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More importantly, to capture the gap of forecasting error from the best modified model, 
the study employs the quantile regression model in line with different conditional 
distributions, say 25%, 50%, 75% and 95%. This guide to analyze the fitted and the 
actual value of the sample toward the residual of forecasting. Since the long term period 
forecasting from 2013m12 to 2017m7 produced the smallest error due to the statistical 
value of RMSE and U index, the study applies this period to predict residuals at 
different conditional quantile. As the result, figure 7 shows that point forecasting 
interval of quantile at 25% and 50% fits perfectly to the actual value and shows the 
errorless rather than the others. 


Table 6: Measurement predictive results from out of sample forecasting 
with n-step ahead 


(1) 0.1386 0.8801 0.1471 0.9562 0.1507 0.9992 0.1470 1.2509 
(2) 0.1322 0.7777 0.1401 0.9498 0.1486 1.1135 0.1547 1.7247 
(3) 0.1338 0.8972 0.1422 1.0534 0.1488 1.0781 0.1464 1.3489 
(4) 0.1394 0.8884 0.1484 0.9416 0.1518 0.9962 0.1482 1.2382 
(5) 0.1325 0.8119 0.1402 0.9458 0.1489 1.1146 0.1555 1.7467 


Source: Author’s estimates 

Note: model | is ARIMA (1, 1, 1), model 2 is ARIMA (3, 1, 4), model 3 is GARCH (1, 1), model 4 is ARIMA (1, 1, 
1) — GARCH (1, 1) and model 5 is an ARIMA (3, 1, 4) - GARCH (1, 1). All models were displayed by time trend 
effect (t). 


Figure 7: Point forecasting interval at different quantile distributions 
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Note: Iny is referred to number of tourist arrivals to Cambodia with logarithm function. 
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Concluding Remarks 


International tourist arrivals are the crucial source of revenue for many developed and 
developing countries. Therefore, to predict its trend is ideal and to manage international 
tourism growth is essential and compulsory to model and forecast adequately tourist 
arrivals and their associated volatility with its order of parameters. This study aims at 
modeling and forecasting tourist arrivals in Cambodia from monthly observations of 
2000m1 to 2017m7, covering 209 samples. The empirical results show that in one part 
toward the baseline regression equation, tourist arrivals have a significant occurrence 
resulted of the appreciation of exchange rate and some internal and external factors such 
as the national election, the AEC, e-Visa application as well as the financial global 
crisis. These factors appearance the power explanation with a statistical significance 
during the observed periods. In another part regarding to autoregressive and volatility 
approach, the empirical results indicate that tourist arrivals is affected by time trend and 
the previous visitors in the past period, say time (t-1). More importantly, the trend is 
reduced due to time lag, say time (t-2, t-3, t-4). The GARCH (1, 1) model of tourist 
arrivals suggests that the short run persistence of shocks lies in the gap of 0.04 while the 
long run persistence lies in the gap of 0.94. RMSE and U index obtained from the 
measurement predictive accuracy of out of sample forecasting reveal that long run 1- 
step ahead of the period 2013m12 to 2017m7 is produced the smallest error among the 
others. Thus, it has more predictive power to apply long term ex-ante forecasting. 


Nonetheless, the empirical findings print out some messages and suggestions for the 
further academic researcher as well as policy makers. Tourism policy makers should 
consider carefully the unexpected events which cause the volatility of tourist arrivals, 
say the economic crisis e.g., and national event such as national election in line with 
some sensitive factors that might disturb strongly to travel decision. In addition, once 
the volatility of tourist arrivals is found, it does matter for tourism and business policy 
makers whether which path of the shock the generating policies to attract tourists could 
be employed effectively. So far, the tourism sector is a relevant economic activity, it is 
significant and necessary to note the unanticipated shock in line with volatility shock, 
will have an inference on tourist arrivals for Cambodia, both in the short and long run. 
Alternatively, it is necessary to determine the extent to which a volatility shock, the 
tourism is diverted to other countries that particularly have similar product development 
of tourism industry. The facts that economic and political events may be changed and 
occurred in the real and exact economic phase with some interaction of an exogenous 
factor. 


References 


Andrew Saayman, Ilse Botha. 2015. "Non-linear models for tourism demand forecasting." 
Tourism Economics 23 (3): 594 - 613. 


Baum, Christopher F. 2015. "ARCH and MGARCH models." EC 823: Applied Econometrics, 
Boston College, Spring 2014. 


Bloomfield, Peter. 2009. "Modeling The Variance of a Time Series." July 31, 2009 / Ben 
Kedem Symposium, Department of Statistics, North Carolina State University. 


Bollershev, T. 1986. "Generalized Autoregressive Conditional Heteroskedasticity ." Journal of 
Econometrics 31: 307-327. 


16 


Chhorn T. and Chaiboonsri C. 


/ Journal of Management, Economics, and Industrial Organization, Vol.2 No.2, 2018, pp.1-19. 


Box and Pierce, D.,A., G. E.,. 1977. "Distribution of Residual Autocorrelation in Autoregressive 
Integrated Moving Average Time Series Models." Journal of the American Statistical 
Association 64. 


Box, G. E. P., and G. M. Jenkins. 1970. "Time Series Analysis, Forecasting and Control." 
Holden - Day. San Francisco, CA. 


Chang, Chia-Lin and Michael McAleer. 2009. "Daily Tourist Arrivals, Exchange Rates and 
Volatility for Korea and Taiwan." Korean Economic Review 25: 241-267. 


Chen, Kuan-Yu. 2011. "Combining linear and nonlinear model in forecasting tourism demand." 
Expert Systems with Applications 38 (8): 10368-10376. 


Chen, Mei-Yuan. 2013. "Time Series Analysis: Conditional Volatility Models." Department of 
Finance, National Chung Hsing University . 


Chia-Lin Chang, Michael McAleer and Dan Slottje. 2009. "Modelling International Tourist 
Arrivals and Volatility: An Application to Taiwan." Department of Applied Economics, 
National Chung Hsing University, Taiwan . 


Chia-Lin Chang, Thanchanok Khamkaew, Roengchai Tansuchat and Michael McAlee. 2011. 
"Interdependence of International Tourism Demand and Volatility in Leading ASEAN 
Destinations." Tourism Economics 17 (3): 481 - 507. 


Chikobvu, Makoni Tendai and Delson. 2017. "Modelling international tourist arrivals and 
volatility to the Victoria Falls Rainforest, Zimbabwe: Application of the GARCH 
family of models." African Journal of Hospitality, Tourism and Leisure 6 (A). 


Chu, Fong-Lin. 2009. "Forecasting tourism demand with ARMA-based methods." Tourism 
Management 30 (5): 740-751. 


Crouch Geoffrey I. 1993. "Currency Exchange Rates and the Demand for International 
Tourism." The Journal of Tourism Studies 4: 45- 53. 


Dickey, D. and W. Fuller. 1981. "Likelihood Ratio Statistics for Autoregressive Time Series 
with a Unit Root." Econometrica (Econometrica) 49: 1057-1072. 


Elzbieta F. and Mirostaw Ga, sowski. 2004. "Modelling Stock Returns with AR - GARCH 
Process." SORT 28 (1 ): 55 - 68. 


Engle, R. F.,. 1982. "Autoregressive Conditional Heteroskedasticity with Estimates of the 
Variance of United Kingdom Inflation." (Econometrcia ) 50 (4 ): 987 - 1008. 


Gang, Jianhua. n.d. "Advanced STATA with Time Series Data 1." Beihang University. 


George Agiomirgianakis, Dimitris Serenis and Nicholas Tsounis. 2014. "Exchange Rate 
Volatility and Tourist Flows into Turkey." Journal of Economic Integration 29 (4): 700 
- 725. 


Ghana | Bennet Otoo. 2016. "Impact upcoming elections on Tourism and Hospitality industry." 
www.myjoyonline.com, November 30. 


Jansky, I., Rippel, M. 2011. "Value at Risk forecasting with the ARMA-GARCH family of 
models in times of increased volatility." IES Working Paper 27/2011 (Charles 
University Marek). 


17 


Chhorn T. and Chaiboonsri C. 


/ Journal of Management, Economics, and Industrial Organization, Vol.2 No.2, 2018, pp.1-19. 


Jewoo Kim, Choong-Ki Lee & James W. Melted. 2016. "Impact of economic policy on 
international tourism demand: the case of Abenomics.”" Current Issues in Tourism 0 (0). 
https://doi.org/10.1080/13683500.2016.1198307. 


José F. P-R. et al., José Francisco Perles-Ribes, Ana Belén Ramon-Rodriguez, Martin Sevilla- 
Jiménez and Antonio Rubia. 2016. "The effects of economic crises on tourism success: 
an integrated model." Tourism Economics 22 (2): 417-447. doi:10.5367/te.2014.0428. 


Kudret Gul, Nuran Aksit Asik, Ali Kemal Gurbuz. 2014. "The effect of global economic crisis 
on Turkish tourism demand and a review for the period 2003-2013." Journal of World 
Economic Research 22-32. doi:10.11648/j.jwer.s.2014030601.14 . 


Kunst, Robert M. 2012. "Econometric Forecasting." Institute for Advanced Studies Vienna and 
University of Vienna . http://homepage.univie.ac.at/robert.kunst/progpres.pdf. 


Louis, Faruk Balli & Rosmy Jean. 2015. "Modelling the tourism receipt’s volatility." Applied 
Economics Letters 22 (2). 


Lukas Falat, Zuzana Stanikova, Maria Durisova, Beata Holkova, Tatiana Potkanova. 2015. 
"Application of Neural Network Models in Modelling Economic Time Series with Non- 
constant Volatility." Business Economics and Management 2015 Conference, 
BEM2015. Procedia Economics and Finance 34 ( 2015 ). 600 — 607. 


Mamula, Maja. 2015. "Modelling and Forecasting International Tourism Demand — Evaluation 
of Forecasting Performance." International Journal of Business Administration 6 (3): 
102 - 112. 


McAleer, Riaz Shareef and Michael. 2005. "Modelling international tourism demand and 
volatility in small island tourism economies." /nternational Journal of Tourism 
Research 7 (6): 313-333. 


MoT. 2012. Tourism Development Strategic Plan 2012 - 2020. Minitry of Tourism. 
http://www.tourismcambodia.org/images/mot/legal_documents/tourism_development_s 
tategic_plan_2012_2020_english.pdf. 


NCSS Ch. 470. n.d. The Box-Jenkins Method. Chapter 470, NCSS Statistical Software, 470-1 
470-14. 


Pathnadabutr, Aliwassa. 2013. "AEC, tourism drive investment trends." Nationa Multi Media, 
February 15. http://www.nationmultimedia.com/business/AEC-tourism-drive- 
investment-trends-30200022.html. 


Paul Leung, Terry Lam & Simon Wong. 2017. "Tourism Development in Cambodia: An 
Analysis of Opportunities and Barriers." Department of Hotel and Tourism 
Management, Hong Kong Polytechnic University, Kowloon, Hong Kong. 


Phillips, P.C.B. and P. Perron. 1988. "Testing for Unit Roots in Time Series Regression." 
Biometrika (Biometrika) 75: 335-346. 


Ramos, Célia M.Q. and Paulo M. M. Rodrigues. 2014. "Tourism Demand Modelling and 
Forecasting: An Overview." Revista de Turismo Contempordaneo 2 (2): 323-340. 


Reider, Rob. 2009 . "Volatility Forecasting I: GARCH Models ." 


Robert R. Andrawis, Amir F. Atiya and HishamEl-Shishiny. 2011. "Combination of long term 
and short term forecasts, with application to tourism demand forecasting." International 
Journal of Forecasting 27 (3): 870-886. 


18 


Chhorn T. and Chaiboonsri C. 


/ Journal of Management, Economics, and Industrial Organization, Vol.2 No.2, 2018, pp.1-19. 


Serenis, Dimitrios and Nicholas Tsounis. 2014. "The Effects of Exchange Rate Volatility on 
Sectoral Exports. Evidence from Sweden, the U.K. and Germany." /nternational 
Journal of Computational Economics and Econometrics. 


Shan, Nada Kulendran & Jordan. 2002. "Forecasting China's Monthly Inbound Travel 
Demand." Journal of Travel & Tourism Marketing 13 (1). 


Shira, Dezan. 2016. "Foreign Exchange Volatility in ASEAN: Assessing Exposure & Managing 
Risk." ASEAN Briefing, March 15. 
https://www.aseanbriefing.com/news/2016/03/15/5230.html. 


Sjo, Bo. 2011. "Estimation and Testing for ARCH and GARCH." Modelling the volatility of the 
Electrolux stock. 


Song et al., Haityan Song, Gang Li and STEPHEN F. WITT AND BAOGANG FET 2010. 
"Tourism demand modelling and forecasting: how should demand be measured?" 
Tourism Economics 16 (1): 63-81. 


Suhejla Hoti, Carmelo Leon and Michael McAleer. n.d. "International Tourism Demand and 
Volatility Models for the Canary Islands." School of Economics and Commerce, 
University of Western Australia . 


Suliman, Ahmed Elsheikh M. Ahmed and Suliman Zakaria. 2011. "Modeling Stock Market 
Volatility using GARCH Models Evidence from Sudan." International Journal of 
Business and Social Science (International Journal of Business and Social Science) 2 
(23 ): 114 - 128. 


Tourism, Ministry of. 2017. Tourism Statistics. Statistics and Tourism Information Department , 
Phnom Penh, Cambodia: Ministry of Tourism. http://www.mot.gov.kh. 


UNWTO. 2013. Economic Crisis, International Tourism Decline and its Impact on the Poor. 
Madrid: World Tourism Organization and International Labour Organization 
(UNWTO). 


Wang, T. 2016. "Forecast of Economic Growth by Time Series and Scenario Planning Method." 
Modern Economy 7 : 212 - 222. 


Webber, Anthony G. 2001. "Exchange Rate Volatility and Cointegration in Tousim Demand." 
Journal of travel Research 39: 398-405. doi:10.1177/004728750103900406. 


Witt, Johann du Preez and Stephen F. 2003. "Univariate versus multivariate time series 
forecasting: an application to international tourism demand." International Journal of 
Forecasting 19 (3): 435-451. 


Yap, Ghialy Choy Lee. 2012. "An Examination of the Effects of Exchange Rates on Australia's 
Inbound Tourism Growth: A Multivariate Conditional Volatility Approach." 
International Journal of Business Studies 111-132. 


Yu-ShanWang. 2009. "The impact of crisis events and macroeconomic activity on Taiwan's 
international inbound tourism demand." Tourism Management 30 (1): 75-82. 
https://doi.org/10.1016/j.tourman.2008.04.010. 


19 


